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Abstract 

t^- ■ It is studied the model on J /if) production allowing to extract parton distribution 

functions (PDFs) from the combined analysis with both data on Drell-Yan and J /if) 
\ production processes. It is shown that this, so attractive from theoretical point of view, 

model, can be safely used in the low energy region E < lOOGeV . The significance of gluon 
contributions to the J /if) cross-sections is investigated. The obtained results in the high 

Q \ energy region occur to be rather surprising. 

\o • 

(N . PACS: 13.65.Ni, 13.60.Hb, 13.88.+e 

^ ■ Nowadays we see the growing interest [H [2] to the close analogy (duality) between Drell-Yan 
^ (DY) H X H 2 -> 7*A -> l+l~X and J /if) E X E 2 -»• J/ipX -»• l+l~X production mechanisms (see 
gn| textbook [3] for review). It is assumed that a such analogy /duality occurs at relatively low 
r-| \ energies, when the gluon-gluon fusion (gg) mechanism of J /if) production is dominated by the 
quark- ant iquark fusion (qq). Then, since J /if) is a vector particle like 7 and the helicity structure 
(N ■ oiqq(J/if)) and (qq)^* couplings is the same, one can get the J /if) production cross-section from 
the DY process cross-section applying the simple replacement 
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where M 2 = Q 2 is the squared mass of dilepton pair, Mj^ ~ 9.59 GeV 2 is the squared J /if) 
O mass and Tj/^ is the full J /if) width. It is believed that the model §B) can be applied in both 
j> \ unpolarized [2] and polarized [1] cases. The later is due to the identical helicity and vector 
^ I structure of 7* and J/ if) elementary channels (all 7 M couplings). In particular, the model given 
^ ' by Eq. <$T§ was exploited in Ref. pQ in the double-polarized (both hadrons in initial state are 
transversely polarized) case. 

It should be noticed here, that though the pure vector 7 M coupling (Dirac coupling) is very 
likely in the J /if) production, still it is a hypothesis. In principle, still the effects due to 
coupling (Pauli coupling) are possible. Later, for the number of processes with the different 
unpolarized colliding hadrons it will be shown that the model (pQ), based on pure vector coupling, 
indeed works well at low energies. It can be considered as a strong argument in favor of pure 
vector coupling mechanism. At the same time, in polarized case we still have no necessary 
data to check the pure vector coupling hypothesis. Thus, for the comprehensive verification of 
the "duality" model it is of crucial importance to collect the data on J /if) production with 
the polarized colliding hadrons. This is a strong additional motivation to perform the such 
experiments. At present the J /if) production processes with polarized proton/deutron target 
are planned at J-PARC [1] (proton beam) and COMPASS [5 J (pion beam). The respective 
experiment with the polarized antiproton beam is possible at GSI (this possibility is now under 
discussion [6]). 

The advantage of model ([T]) is that in the region of w-quark dominance (large Bjorken x) 
all couplings exactly cancel out in the ratios of cross-sections (like asymmetries), so that they 
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become absolutely the same for the DY and J /if) production processes (see Eq. (23) in Ref. [T] 
and Eq. (10) in Ref. [2]). Thus, in this kinematical region it does not matter for the u quark 
PDFs extraction where do the dilepton pair production events come from: from continuum 
or from J /if) production region. Certainly, the such possibility to use J /if) production for 
PDFs extraction is very attractive because the dilepton production rate in the J /if) production 
region is two orders of magnitude higher than in the continuum region above the J /if) mass. 
In particular, the model ([I]) give us the possibility to extract transversity h\ u and the first 
moment of Boer-Mulders function hf^ from J /if) production region. Namely, in the large x 
region the equations (Eqs. (19), (20) in Ref. [7J and Eqs. (12), (13) in Ref. [8]) for unpolarized 

#i# 2 -> 7*x -> i+rx 

flu{xA)\Aflu{x B )\B + (U <-> U) 

and sing le-polarized H^j -> 7 *X -»■ l+l~X 

£ _ 1 h m l) ( x a) \Ah lu (x B i ) \bT_ + (u <-> u) 
2 fm(x A )\ A f lu (x B ^)\ B T + (u ^ u) 

DY processes are the same for the respective J /if) production processes HiH 2 — > J/if)X — > 
l + l~X and HiH\ — > J/if>X — > Z + Z~X, providing us an access to and as a result of 
combined analysis of DY and J / if) data. 

Certainly, the "duality" model ([1]) is applicable only in the such kinematical regions where 
among the elementary processes contributing to J /if) production the quark-antiquark fusion 
process dominates over the gluon-gluon fusion. Qualitatively, it is clear that the gluon contri- 
butions should be suppressed at low energies: at fixed Q 2 value (9.6 GeV 2 for J /if) production 
we deal with) x increase when s = (p Bl + Ph 2 ) 2 decreases, while the ratio g(x)/q(x) sharply 
decreases when x increases. However, we certainly need the quantitative analysis to find the 
region where we can safely neglect the gluons. The main goal of this paper is to perform the 
such quantitative analysis. To this end we, besides of the model (jTJ), will consider the most 
popular and well established model on J / if) production which includes all elementary processes 
(quark-antiquark as well as gluon-gluon fusion) and believed to be valid at any energies of 
incident hadrons - see [9] and references therein. Then we will compare the predictions of both 
these models with the respective existing experimental data |10j-|16j. 

To cancel unknown constants, we will study not absolute cross-sections but the ratios^ of 
cross-sections corresponding to the variety of hadrons/nucleus in initial state. Namely, we will 
consider the ratios of the angle and xf integrated (the later in the forward hemisphere xp > 0) 
cross-sections on J /if) production <J pp /cr n ±p, (Tpa/p-k^a and a pp /ap p , where the symbol A denotes 
the different target nucleus (these are W, C, Ca, Cu, Pt here). The point is that the hadrons 
7i ± and p (on the contrary to proton and nuclei) contain the antiquark in the valence state. 
That is why it is very important to study the ratios like <J pp /cr n ± p and c pp /a pp which should 
possess very specific behavior: they should sharply decrease when s decrease (large x, valence 
quark/ antiquark dominance region) and they should increase up to unity when s takes large 
values (small x, sea quark and gluon dominance region). The most of data on J /if) production 
giving us an access to the such kind of ratios are the data in the forward hemisphere with the 
pion or antiproton beams colliding with the proton or nuclear targets |10j-[16j. 



1 Let us recall that as a final goal we are interested namely in the ratios - asymmetries, which can give us 
an access to the different PDFs. 
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The simple "duality" model ([T]) applied to the ratio of cross-sections [(Th^/ch' H 2 ]x F >o 
gives 

1 — 771^ Is 

r , , Jo «*' dx F [s(x 1 +X 2 )]- 1 F**> 

where the quark/antiquark flow F^ lH2 is given by 

F S lH2 = E [? Hl (^)ff ffa (^)+ff ffl (^i)5 ffa (^)], (5) 

q=u,d,s 



and Xi j2 are expressed via x F = X\ — x 2 as x li2 = [±x.f + yjxp + 4m 2 ^/ s]/2. Notice that 

obtaining Eqs. (j4]), (JSJ) we put the ratios g J d ' / ' g J J^ and g J J^ ' j ' g J J^ equal to unity in accordance 
with the existing experimental data [T7]. Indeed, the data [TF] on J/ip production with the 
absolutely symmetric (so that the cross section per nucleon reads cr n ±c = [0V±p + 0"7r±n]/2) 
carbon target clearly indicates that the ratio <J n +c / 'o'n-c is seem to be near unity in J /if) 
region, while it falls toward 1/4 above the dilepton pair mass - see Fig. 2 in Ref. [T7] and 
discussion around. This is a good argument that on the contrary to qq annihilation mechanism 
for DY process (where d quark is suppressed by the factor e 2 Je 2 u = 1/4), for J ftp production 
u and d quarks should enter the cross-sections with the same charge factors: g J J^ '/ ' g J J^ — 1. 
We also use the analogous relation g J J^ ' / ' g J J^ — 1 keeping in mind that the squared sea strange 
quark contribution is not so significant. 

The main difference of the "gluon evaporation" model from the model ([T]) is that the former 
in addition to qq fusion contribution F^ lH2 given by Eq. (jSJ) contains also gluon-gluon fusion 
contribution F^ lH2 = g Hl (xi)g H2 (x2). Applied to the ratios [ch^I (J h' 1 h 2 ]x f >o-i the "gluon 
evaporation" model produces: 



O~H 1 H 2 \x F >0 _ {(?qq + (^gg) H 1 H 2 \x F >0 
O~H[H 2 \x F >0 (Cqq + O'gg) H[H 2 \x F >0 



(6) 



HiH 2 
99(99) 



Am 2 rl-— 



. °dQ 2 s dx F a^ c5 ^ c5 \Q 2 ) „ 9/ XlX2 F^, (7) 

*f>0 J 4m 2 V J * ^ J Q 2 (x 1 + X 2 ) qq (99)^ V7 



where 2m c = 3.0 GeV and Imp, = 3.7AGeV are respectively the cc and open charm thresholds, 
while the elementary cross-sections a qq ~* cc , a"^ cc are proportional to a s (Q 2 ) and can be found, 
for example, in Ref. [9] (see Eqs. (3) and (4)). For the comparison purposes we will consider also 
the "gluon evaporation" model without gluon contribution F^ g lH2 . It is obvious that it differs 
from the "duality" model only by the extra integration over Q 2 with the weight a qq ^ cc (Q 2 ) . 

Let us first consider the ratios app/a n ± p . The results obtained within the "duality" and 
"gluon evaporation" models are presented in Fig.l in comparison with experimental data. First 
of all, one can conclude that in the low energy region, near the first experimental point y/s ~ 
8.7 GeV, the curves corresponding to "duality" model and to "gluon evaporation" model with 
and without gluons almost merge and equally well describe the existing experimental data. 
This is not suprising and is in agreement with the qualitative predictions: as it was discussed 
above, the gluon contribution should be suppressed in the low energy region. At the same time, 
the results in the high energy region occur to be rather surprising: even at very high energies 
(150 GeV and 200 GeV) the gluon contribution seems to be insignificant in the ratios (Tpp/cr^p 
and the curves with and without gluon contributions F^ lH2 equally well describe the existing 
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data. Thus, in the case of 150-=-200 GeV pion beam and the proton target we need to improve 
the quality of data on J/ijj production to distinguish between the models with and without 
gluon contribution. The such measurements can be performed, for example, by COMPASS 
collaboration [5], where the possibility to study Drell-Yan and J /if) processes with the pion 
beam is now under consideration. 

The absolutely analogous picture, insignificance of the gluon contribution even at high 
energies, occur also for the ratios v v aI&-k±a with the different target nuclei - see Fig. 2. In Fig. 2 
the data with the approximately equal ratio Z/A ~ 0.4 are collected. The curves corresponding 
to the model calculations are obtained with Z/A = 0.4 and neglecting^] the nuclear effects, so 
that the cross-section per nucleon reads: (ThA\h=-K±,p — \ a h v + U — -f J &hn- 

Let us now consider the ratios o vv /a pv with the incident antiproton instead of pion. The 
results are presented in Fig.3. While in the low energy region we again see the good agree- 
ment between the models with and without gluons and the data (as it should be from the 
qualitative consideration), the situation in high energy region is absolutely different. First, one 
can conclude, that the gluon contribution becomes very significant in this kinematical region. 
Second, and rather surprising conclusion, is that the widely used "gluon evaporation" model 
works rather bad in this case - the respective curve lies well below the data (bold solid line in 
Fig. 3). Notice that this result is in a strong disagreement with the statement made in Ref. 
[15J, where the same experimental points were used. Indeed, calculations within the "gluon 
evaporation" model presented in Ref. [15] (bold dashed line in Fig.3) are in a good agreement 
with the data points which was claimed there as a strong argument in favor of that model (see 
Fig. 6 in Ref. [15] and discussion around). The reason of this disagreement is in the gluonEl 
sector of the model, because the main difference of the parametrization [18] used in Ref. [15] 
and the modern parametrizations^] used in our calculations is just the gluon distribution values. 
Certainly, one should trust namely the calculation where the modern parametrizations are used, 
because G(x) is much better fixed there by a lot of new data appeared since the paper [18] issue. 
The respective results occur to be in the strong disagreement with the data and, thus, one can 
conclude that to pass this test the gluon sector in the "gluon evaporation" model should be 
essentially modified. Besides, all existing nowadays models on J / tjj production should also pass 
this test on high energy behavior, and this is a subject of our future investigation. 

In conclusion, we have tested the "duality" model as well as "gluon evaporation" model in 
the different energy ranges. It is shown that the "duality" model (as well as "gluon evaporation" 
model) works well in the low energy sector, s < lOOGeV 2 . In this region the curves with and 
without gluon contributions almost coincide and equally well describe the data. Thus, we can 
safely apply in this kinematical region so attractive from the theoretical point of view "duality" 
model. This give us the unique possibility to use the J/tp production region for combined 
analysis with the data on Drell-Yan processes to extract PDFs we are interested in. This is an 
excellent possibility to essentially increase the statistics since the J/ ip production cross-sections 
are in about two orders of magnitude higher than the respective DY cross-sections. 

On the other hand, we got two surprises in the high energy region. First is for the pion- 

2 Usually [§] the nuclear effects are accumulated in the multiplier A a . However, for xp integrated cross- 
sections these factors differ a little from unity - see Ref. [9] and references therein. 

3 It is clear seen from Fig. 3. Indeed, the curves corresponding to the "gluon evaporation" model without 
gluons, for both old and modern parametrizations (thin dashed and solid lines in Fig. 3), practically merge. 

4 We present here the result (bold solid line in Fig. 3) with the modern and widely used GRV98 [TH] 
parametrization. However, our calculations with another modern parametrizations produce the same picture 
(the results differ a little from the respective results with the GRV98 parametrization). 
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proton/nuclei collisions, where the gluon contribution seems to be insignificant even at the 
energies of incident pion about 150-^200 GeV. The second surprise occurs for the antiproton- 
proton collisions. Here, on the one hand, the gluon contribution in the ratio <J vv /a pv is very- 
essential. On the other hand, it is not properly described by the popular and widely used "gluon 
evaporation" model, when we use the modern parametrizations on G(x) instead of old (and 
rather incorrect) one. Thus, it seems that the "gluon evaporation" model should be properly 
modified. 

Thus, the results of performed tests definitely show that both theoretical and experimental 
efforts are still necessary to answer the raised questions. We need the new information from 
both low and high energy regions. In this respect it should be noticed that at present the new 
experiments with the proton-proton collisions on Drell-Yan and J /if} processes are planned [lj 
at J-PARC facility. The respective measurements with the low energy antiproton beam could 
be performed at GSI (see Ref. [6]). At the same time, the researches on Drell-Yan and J /if) 
physics at high energies (150-^200 GeV energy of incident pion) are now planned at COMPASS 
experiment [5]. 
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a(pp -> J/i)/)/o(7t'p -> J/y) (pp _^ J/ v )/ ( re * p _> J/ v ) 




Figure 1: Ratios of cross-sections cr pp /a^ p on J/ip production calculated with two models in 
comparison with the experimental data. Solid line corresponds to the "duality" model, Eqs. 
(HI), (EJ). Dashed line corresponds to the "gluon evaporation" model, Eqs. (J6]) , (JTJ) . Dot-dashed 
line corresponds to "gluon evaporation" model without gluon contribution. Here GRV94 [21] 
parametrization for the proton PDFs and GRV [20] parametrization for the pion PDFs are 
used. The experimental data (points with error bars) are taken from [9] (Tables 2 and 3). 
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Figure 2: Ratios of cross-sections o- p a/o-^a on J/i> production for the nuclear targets 
(Z/A ~ 0.4) calculated with two models in comparison with the experimental data. Solid 
line corresponds to the "duality" model, Eqs. ([4]), ([5]). Dashed line corresponds to the "gluon 
evaporation model", Eqs. ([6]), ([7]). Dot-dashed line corresponds to the "gluon evaporation" 
model without gluon contributions. Here GRV94 [21] parametrization for proton PDFs and 
GRV [20] parametrization for pion PDFs are used. The points with error bars correspond to 
experimental data. First point: W, Z/A=0.40 [IT] ; second and third points: Pt, Z/A=0.40 
P3|; fourth point: C, Z/A=0.5 |J6]; fifth point: Be, Z/A=0.44 HT 
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Figure 3: Ratios of cross-sections a pp /ap p on J /if) production for the proton target calculated 
with two models in comparison with experimental data. The experimental data are taken from 
Ref. [15], where the respective data of NA3, WA39 and UA6 collaboations were collected (see 
Fig. 6 in Ref. P3]). The bold solid and dashed lines correspond to calculation with the "gluon 
evaporation" model (with GRV98 [19] and Duke-Owens [H] parametrizations, respectively). 
The thin solid and dashed lines correspond to calculations with the "gluon evaporation" model 
without gluon contrubutions (with GRV98 [19] and Duke-Owens [18] parametrizations, respec- 
tively). The dotted line correspond to calculation with the "duality" model and GRV98 [T9] 
parametrization. 
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